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Due to their technological importance, point defects in silicon are among the best studied physical 
systems. The experimental examination of point defects buried in bulk is difficult and evidence for 
the various defects usually indirect. Simulations of defects in silicon have been performed at various 
levels of sophistication ranging from fast force fields to accurate density functional calculations. 
The generally accepted viewpoint from all these studies is that vacancies and self interstitials are 
the basic point defects in silicon. We challenge this point of view by presenting density functional 
calculations that show that there is a new fourfold coordinated point defect in silicon that is lower 
in energy. 

PACS number: 61.72.Ji 



The stability of crystalline silicon comes from the fact 
that each silicon atom can accommodate its four valence 
electrons in four covalent bonds with its four neighbors. 
The traditional point defects in silicon, the vacancy and 
the various interstitials, are obtained by taking out or 
adding atoms to the crystal and thus all destroy fourfold 
coordination. A relatively high defect formation energy 
for these defects is the consequence. In addition there 
is a point defect, that conserves the number of particles, 
the Frenkel pair, consisting of a vacancy and an inter- 
stitial. If the vacancy and the interstitial are close, the 
formation energy of a Frenkel pair is less than the sum of 
an isolated vacancy and interstitial (Table |). Neverthe- 
less the formation energy is still considerable since again 
bonds are broken. Figure |l| shows a novel defect config- 
uration that has, in contrast to all other point defects, 
perfect fourfold coordination and will therefore be called 
Four Fold Coordinated Defect (FFCD). 

Using state of the art plane wave density functional 
theory (DFT) calculations we will now present evidence 
that the formation energy of the FFCD of 2.4 eV is lower 
than the formation energy of all other known point de- 
fects both in intrinsic and doped silicon. Even though 
several calculations of this type were published for the 
traditional defects, we have decided to repeat them for 
several reasons, a) All these calculations take advan- 
tage of error cancellations to obtain energy differences 
that are more precise than the total energies themselves. 
This cancellation is obviously best if all the total ener- 
gies that are compared are calculated with exactly the 
same method, b) Most DFT calculations used the basic 
Local Density Approximation (LDA) whereas we used a 
more precise General Gradient Approximations 0]. c) 
As was recently shown ^ cells of at least 216 atoms 
are required for a reasonable convergence. Most pub- 
lished calculations were done with smaller cells, d) We 
constructed very accurate pseudopotentials ^ based on 
atomic calculations with the same density functional as 



used in our target calculation and we afforded a very 
large plane wave basis set (25 Ry for the LDA and 35 Ry 
for the GGA calculations) . e) The point defect formation 
energies in doped silicon were up to now obtained in an 
approximate way by combining density functional total 
energy results with the concept of a chemical potential 
that varies as a function of doping across the experimen- 
tal band gap of 1.12 eV of silicon |||,^. We have modeled 
the effect of doping in a coherent way for both n and 
p type doped silicon by including explicitly two doping 
atoms (aluminum and phosphor) in our 216 atom su- 
percell. This procedure leads of course to high doping 
concentrations for any reasonable size of our computa- 
tional cell, but by treating coherently and accurately the 
extreme cases of high and zero doping we can expect the 
numbers for more moderate doping concentrations to lie 
in between these results. 

Table | shows all the total energy results. We have in- 
cluded in our compilation all defects that have been clas- 
sified as being low in energy in previous calculations [Q-Q , 
i.e. the split (110) X interstitial, the hexagonal H inter- 
stitial and the vacancy. In addition we have included a 
low energy Frenkel pair. In agreement with recent calcu- 
lation by Needs , we found that both the tetrahedral 
interstitial and the caged interstitial are higher in 
energy and meta-stable. They are therefore not included 
in the compilation. The energies in Table || do not con- 
tain an additional lowering of the energy that is obtained 
in an infinite crystal by a long range elastic field, which 
is suppressed in our finite cell. We have calculated this 
extra elastic energy for a few configurations and found 
it to be negligible, namely less than .05 eV. Because of 
the negative U character of the traditional vacancies and 
interstitials that was postulated with the help of DFT 
calculations 0,|2j and confirmed experimentally [p^ , 
the lowest energies are obtained by closed shell config- 
urations. Consequently all calculations were performed 
without spin polarization. 
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For the FFCD the bond length and angle do not signif- 
icantly deviate from their bulk values. Whereas the bond 
length and angle in the bulk are 2.35 A and 109 degrees, 
they respectively vary from 2.25 to 2.47 A and from 97 
to 116 degrees for the bonds formed by the two (red) 
defect atoms. The fourfold coordination is also visible 
from the fact that the centers of the maximally localized 
Wannier functions |jl^ , that give the best possible spatial 
location of an electron, are located practically exactly in 
the middle of the bonds established by the geometric dis- 
tance criteria. The spread of all the Wannier functions 
is very close to the bulk value of 2.7 a.u. The fact that 
the LDA and GGA values are very close is not surprising. 
The simple covalent bonds are very well described by any 
reasonable density functional and we expect an error of 
less than a tenth of an eV for this configuration. 

Correlation effects are more important for the broken 
bonds of the other point defects. Hence the DFT values 
are less reliable than in the case of the FFCD. Recent 
Quantum Monte Carlo(QMC) calculations |^ for inter- 
stitials are in much better agreement with GGA values 
than with LDA values. We therefore believe that it is 
unlikely that the energetical ordering predicted by our 
GGA calculations does not correspond to reality. 

The broken bonds of the other point defects are again 
best visualized by looking at the centers of the maximally 
localized Wannier functions in the neighborhood of de- 
fects. They are located in empty space regions rather 
than at positions where geometric criteria would pre- 
dict bonds. Their spread is also significantly larger than 
the bulk value. In the case of the vacancy the two non- 
bonding Wannier functions have a spread of 3.6. a.u., for 
the X interstitial 4.1 a.u. and for the H interstitial 4.7 
a.u.. For the Frenkel pair there are two non-symmetric 
Wannier functions with spreads of 3.5 and 4.4 a.u.. 

Whereas the classical point defects were all postulated 
based on simple symmetry considerations, the discov- 
ery of this new point defect was made possible by new 
algorithmic developments. The configurational space 
was systematically explored by using a modified basin 
hopping method and an inter-atomic silicon poten- 
tial (l|]. The most promising configurations were then 
refined with a plane wave electronic structure code p^ ] 
developed by J. Hutter et al.. The preconditioned con- 
jugate gradient method ||l^ was used for the electronic 
optimization and the DIIS method | |2C| ] without any sym- 
metry constraints for the ionic relaxation. 

The remarkable fact that the FFCD has not yet been 
detected experimentally, in spite of its equilibrium con- 
centration which should be many orders of magnitude 
larger than that of the other defects, is probably due 
to two factors. First, the experimental search for de- 
fects was always guided by theoretical predictions. Sec- 
ond, it is invisible with standard experimental tech- 
niques. Because of the perfect fourfold coordination no 
unpaired electrons exist, a necessary condition for Elec- 



tronic Paramagnetic Resonance (EPR) and Electron Nu- 
clear DOuble Resonance (ENDOR) experiments. Deep 
Level Transient Spectroscopy (DLTS) requires electronic 
levels within the gap. Given the perfect coordination 
this seems unlikely. To examine this point further we 
have calculated the LDA Highest Occupied Molecular 
Orbital (HOMO), Lowest Unoccupied Molecular Orbital 
(LUMO) splitting. Even though it is known that the 
gap is not well reproduced by DFT, these levels allow 
for a qualitative interpretation. The results (Figure H) 
show that the FFCD disturbs the band structure indeed 
much less than the other defects. The FFCD spectrum 
is also insensitive to the doping level whereas the other 
spectra are strongly influenced via various Jahn Teller 
distortions ||l^ because of their orbital degeneracy. 

The basic FFCD can serve as a building block for more 
extended defects. Two simple defects obtained by com- 
bining just two FFCD's in two different ways are shown 
at the bottom of Figure I. Their energies of 3.98 eV for 
the tilted configuration and of 4.17 eV for the parallel 
configuration are significantly lower than the sum of the 
energies of two isolated FFCD's. 

The consequences of our results are wide ranging. A 
reexamination of numerous experimental results will be 
necessary. In particular the tacit assumption that the 
same kind of defects are responsible both for the elec- 
trical and diffusion properties of silicon has been refuted 
by our calculations. Even though the FFCD does pre- 
sumably not strongly influence transport properties, it is 
expected to play an important role in the diffusion prop- 
erties of silicon. 

We acknowledge the most useful comments of Giovanni 
Bachelet, Erik Koch, Richard Needs and Cyrus Umrigar 
on our manuscript. 
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TABLE I. The GGA formation energy in eV of the point defects described in the text and depicted in Figures |l|and |2|For 
comparison the LDA values are given in parentheses for intrinsic silicon. 





p-type 


intrinsic 


n-type 


FFCD 


2.45 


2.42 (2.34) 


2.39 


Frenkel 


5.65 


4.32 (4.26) 


5.77 


X interstitial 


3.33 


3.31 (2.88) 


2.98 


H interstitial 


2.80 


3.31 (2.87) 


3.12 


Vacancy 


3.01 


3.17 (3.56) 


3.14 
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FFCD defect 



FFCD defect 



Two non-parallel FFCD's Two parallel FFCD's 

FIG. 1. Two views of the FFCD defect (upper part) and two FFCD pairs (lower part). The defects are obtained by moving 
atoms from the initial positions denoted by black spheres into the final red positions. The centers of the Warmier functions are 
indicated by the small blue spheres for the single FFCD but omitted for the two pairs. The formation of the FFCD is easily 
understandable from the upper left panel. There, the upper red atom originates form the upper black position to the left. The 
swinging movement needed for this displacement does not break the bonds with the upper neighbors, but with the lower ones. 
A corresponding swing is needed for the lower defect atom. In the final position two new bonds can be formed so that in the 
end all atoms are again fourfold coordinated. 
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X interstitial 



H interstitial 



Vacancy Prenkel pair 

FIG. 2. Classical point defects in silicon: Red atoms indicate atoms in new positions. For the vacancy and the Frenkel pair 

the empty lattice site is shown by a black sphere. The centers of extended Wannier functions are visualized by increasingly 
bigger blue spheres, Wannier functions with the typical bulk extension by small blue spheres. 
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BULK FFCD X Intst. H Intst. Vac. Frenkel 



FIG. 3. The HOMO LUMO splitting for the configurations of Table |. The HOMO is denoted by a solid fine, the LUMO by 
a dashed line. Red corresponds to intrinsic silicon, green to p-doped and blue to n-doped silicon. The two levels of intrinsic 
silicon without defects are drawn over the whole interval. 
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This figure "INT_h.jpg" is available in "jpg" format from: 
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